Adiabatic Quantum Computing in systems with constant inter-qubit couplings 
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We propose an approach suitable for solving NP-complete problems via adiabatic quantum com- 
putation with an architecture based on a lattice of interacting spins (qubits) driven by locally 
adjustable effective magnetic fields. Interactions between qubits are assumed constant and instance- 
independent, programming is done only by changing local magnetic fields. Implementations using 
qubits coupled by magnetic-, electric-dipole and exchange interactions are discussed. 
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I. INTRODUCTION 

It is one of the central open questions in the field 
of quantum computing if quantum algorithms can effi- 
ciently solve computationally hard problems in combi- 
natorial optimization. The most common problems of 
this kind encountered in practice belong to a so-called 
NP-complete class p| and are in almost one-to-one cor- 
respondence with the spin glass models in physics |2( . A 
basic property of these models is an onset of an exponen- 
tially large number of deep local minima of the energy 
landscape formed by spin configurations separated by a 
large number of spin flips. 

A general framework for solving optimization prob- 
lems on a quantum computer is provided by an Adiabatic 
Quantum Computation (AQC) [3j,|J]. In its simplest form 
AQC corresponds to a "quantum annealing" (QA) of a 
spin system in the uniform transverse magnetic field de- 
scribed by a smoothly-varying in time Hamiltonian H 
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H(T) = -rJ2*i + Ho, T = T(t). 



(1) 



r is proportional to the transverse field strength; af is 
a Pauli matrix for the i th spin. The first term in Q is 
called a "driver" , it causes transitions between the eigen- 
states of the problem Hamiltonian Hq whose ground state 
encodes the solution of the classical optimization prob- 
lem in question. For example, to find the ground state of 
the classical Ising model one uses 
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where hi are longitudinal fields acting on spins and the 
summation is over the pairs (i, k) of coupled lattice sites. 
In QA the value of Y = T(t) > is slowly decreasing in 
time from a large value Yo 3> max(|Jy|) at the start, to 
r = at the end of the algorithm. The initial state is pre- 
pared to be a ground state of H(Y ) ~ — ro YljLi &j with 
each spin pointing in the positive x-direction. For adia- 
batically slow variation of Y(t) the system state will be 
closely tracking the instantaneous ground state \tpo(T)} 
of H(Y). It will approach the ground state of Hq at the 



end of the algorithm and the solution of the optimization 
problem can be recovered by measurements. 

A metric for the performance of the AQC can be given 
in terms of the minimum value g m i n = minr.g(r) of the 
excitation gap g(Y) between the ground and first excited 
energy levels of H(Y). In particular, if l/g m i n grows no 
faster then polynomially in the problem size N then so 
does the runtime of AQC. 

It was demonstrated recently for AQC applied to com- 
putationally hard random instances of NP-complete Sat- 
isfiability problem that the algorithm performance is sub- 
stantially affected by the existence of the first-order quan- 
tum phase transition for some value of Y = Y* in the limit 
N — > oo 5] . In this limit the instantaneous gap averaged 
over the ensemble of random problem instances vanishes, 
(g(Y*)) — > 0. This produces a significant difficulty in the 
analysis of the asymptotic complexity of AQC. Indeed, 
assume that for a given problem instance 2jv of a finite 
size N the minimum gap is achieved at some value of 
r = Yj N . For large N the distribution of the values of 
Txjv over the ensemble of random problem instances In 
is peaked around T* and has a width ~ N^ 1 / 2 . There- 
fore the minimum of the ensemble-averaged excitation 
gap, minr(<?), will scale down polynomially with N. At 
the same time the ensemble average of the minimum gap 
(minr g(Y)} can be exponentially small in N. 

The above discussion implies that, in general, the or- 
der of time-minimization (over Y) and statistical averag- 
ing of the excitation gap over the ensemble of problem 
instances cannot be inverted. Therefore static properties 
of the quantum phase transition (e.g., a phase diagram) 
are not sufficient to obtain the true asymptotic complex- 
ity of AQC and a daunting theoretical task of studying 
the dynamics of the phase transition is required. An 
crucial insight here can be gained from the experiments 
implementing AQC for problem instances of a large size. 

An early work on this subject is a quantum anneal- 
ing experiment 0. It uses the macroscopic samples of 
a disordered magnet LiHo x Yi_ x F4 with the dipolar cou- 
pling between the spins that are formed by the doublet 
states of Ho +3 ions randomly substituted for nonmag- 
netic Y 3+ . When the magnetic field is applied perpen- 
dicular to the Ising axes the system can be approximated 
by the Hamiltonian H |JTJ,J2J corresponding to a quan- 
tum 3D Ising model with the random antiferromagnetic 
interaction between neighboring spins in the transverse 
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magnetic field T = T(t). 

For r = finding a ground state of H can be mapped 
onto an NP-complete problem 0. However this map- 
ping is only approximate for the random magnet 6] be- 
cause it is not exactly of the Ising type due to the very 
small components of the g-factor transverse to the Ising 
axis. Also the analog computation via QA [(| if of a 
low-fidelity type because it does not have any explicit 
control nor correction of the decoherence-induced errors, 
unlike the conventional quantum computation schemes 
based on quantum circuits. Finally, the adiabatic evo- 
lution in QA [fj collapses beyond the point r* > of 
a quantum phase transition from a quantum paramag- 
net to a quantum spin glass where the system escapes 
from the adiabatic ground state to states with low-lying 
energy levels forming a quasi-continuous spectrum |8( ■ 

Despite all these shortcomings the magnetic suscep- 
tibility study in ;6j and the subsequent numerical sim- 
ulations of this system Q shows that QA provides a 
dramatic speedup of convergence to low-energy states of 
LiHo x Yi_ x F4 over a purely thermal annealing procedure 
at zero transverse magnetic field. This was attributed 
to the fact that the tunnelling processes between local 
minima of the energy landscape open up new transition 
pathways as compared to thermally activated spin flips. 

The decoherence rate for the lowest doublet (spin) 
states in Ho ions comes from the coupling to Ho nuclei 
and is small [11,^3 (a), possibly much smaller then the 
level splitting and magnetic dipolar coupling between the 
neighboring Ho spins One can expect that the suc- 
cess of QA is based substantially on a strong build-in 
quantum coherence with extremely entangled many-spin 
states and quantum correlations extending out over large 
group of spins (cf. E3(b)). 

All NP-complete problems can be mapped onto one 
another by a classical algorithms that scale polynomi- 
ally in the problem size N lj. However different prob- 
lems correspond to spin Hamiltonians with quite different 
physical properties, such as long-range vs short-range in- 
teractions, fully-connected vs sparsely-connected graphs, 
etc. Therefore it would be of a fundamental importance 
to experimentally study the properties of QA in meso- 
scopic systems modelling different NP-complete problems 
at r = 0. This can shed a light on the power of QA even 
if the experiments are of a low-fidelity type @ . However 
one has to overcome several obstacles to implement this 
approach: 

(i) many NP-complete problems correspond to a long- 
range interaction between classical spins while the 
underlying physical systems have short-range inter- 
actions; 

(ii) in most of the perspective material systems the cou- 
pling between qubits is either difficult or impossible 
to control. The former includes systems with super- 
conducting qubits ElHEllliim and the latter 
includes systems with magneto-dipole [53, [3, HH > 
electro-dipole [l7lll^ . lT^.l20| . and and elastic-dipolc 



0,|22 interactions; 

(iii) irregularities of the crystal lattice and random spa- 
tial positions of qubits give rise to the fluctuations 
in magnitude and sign of the qubit coupling coeffi- 
cients across the sample. 

Kaminsky, Lloyd and Orlando (KLO) [2(| proposed a 
novel AQC architecture for solving NP-complete Max- 
imum Independent Set problem using superconducting 
qubits. In their method an instance of the problem is 
converted into an equivalent circuit with ferromagnetic 
and antiferromagnetic couplings in a uniform magnetic 
field. In |26( switchable interqubit couplings using one 
of the mechanisms presented in |l3l |28| are necessary for 
encoding an instance of MIS, but the architecture does 
not require a change in couplings during the algorithm 
execution for a given problem instance. 

At the same time the couplings in [2g have to be read- 
justed for solving a new instance of MIS problem. Also 
the scalability of this is limited because it only allows 
the switching between the regimes of small and large 
couplings. I.e., it cannot guarantee that two sites are 
completely decoupled, nor the sign of the coupling can 
be changed. 

In what following we propose the architecture for AQC 
that assumes that qubit couplings are not adjustable 
and addresses all the problems (i)-(iii) mentioned above. 
Moreover the architecture that we propose for solving 
various NP-complete problems via AQC is universal, in 
a sense that it is independent on any particular qubit 
implementation. 

II. MAXIMUM INDEPENDENT SET PROBLEM 

For practical implementation of AQC for solving an op- 
timization problem it is necessary to translate this prob- 
lem into a model Hamiltonian with single-spin and pair- 
wise interactions only. It is especially straightforward to 
do for the NP-complete graph theory problem Maximum 
Independent Set (MIS) J2(j. This problem is defined on 
a graph G = (V, E) that is composed of the set of ver- 
tices V connected by edges from the set E. MIS is the 
problem of finding the largest subset S of vertices V such 
that no two vertices on the subset share an edge from E. 
The reason why this problem is well suited form AQC is 
that its solution is isomorphic to finding a ground state 
of an antiferromagnetrically coupled Ising spin model in 
a uniform magnetic field \J\. 

The next obstacle has to do with the fact that we want 
to use spins on a 2D lattice and use only nearest-neighbor 
interactions. If the graph is non-planar it cannot be 
embedded into 2D lattice with nearest-neighbor inter- 
actions only. And there is no guarantee that a graph 
just constructed has a drawing on a plane without inter- 
secting edges. Fortunately using a construction invented 
by Carey, Johnson and Stockmeyer |27| we can replace 
two intersecting edges by a gadget of Fig.^ Maximizing 
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the independent set in the gadget reveals opposite ver- 
tices cannot both be occupied without incurring energy 
penalty. Solving MIS problem for the transformed graph 
gives the solution for the original one. Replacing all in- 
stances of intersecting edges by such gadget renders the 
graph planar. For this reason, without losing generality, 
we can restrict our discussion to planar graphs. 




FIG. 1: Every intersection of two edges is replaced by a planar 
graph (right). 

Before we describe the embedding of planar graph into 
a regular 2D lattice let us demonstrate how the MIS 
problem maps onto Ising model. We model vertices of 
the graph with classical Ising spins that can take values 
Si = ±1. The edges correspond to antiferromagnetic in- 
teractions between spins. Consider the following classical 
Hamiltonian: 
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where J^k < corresponding to antiferromagnetic inter- 
action and local magnetic fields are hi = ^ fc J^fc (the 
sum is over the neighbors of a vertex i) . This expression 
can be (up to a constant) rewritten as 
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Since all interactions are antiferromagnetic (J{ t k < 0) 
and (1 — Si)(l — Sfc) > 0, we must have Hq > 0. The 
minimum value Hq = is reached only if no edge (i, k) 
has both s,; = — 1 and Sk = — 1. In other words the set of 
vertices {i|sj = —1} is independent, no two vertices from 
this set share a vertex. Hence, the minima of Hq are 
all possible independent sets. Next, we adjust individual 
magnetic fields h' i = hi — J, which is equivalent to adding 
a perturbation 
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Its role is to maximize the number of vertices with 

I . i.e. maximize the cardinality of MIS. By choosing 
J < minuty | Jj,*; | we guarantee that all minima of Hq + V 
correspond to independent sets with Hq = 0. Indeed, if 
for some minimum of Hq + V and for some i, k both s, 

1 and Sfc = —1, then setting Sj = +1 decreases Hq by at 
least 4| Ji k | > 4 J and at the same time increases V by 2 J. 



The overall energy decreases by at least 2 J contradicting 
our assumption that we started from a minimum of Hq + 
V. Therefore, the ground states of Hq + V correspond 
to solutions of MIS problem. Moreover, the gap between 
ground state(s) and first excited state(s) is at least 2J. 

Obviously, very few planar graphs can be drawn on a 
lattice, with vertices corresponding to the nodes of the 
lattice and edges corresponding to nearest-neighbor links. 
In general case vertices have to be represented by whole 
clusters of spins. We generalize Hamiltonian |J5J| as fol- 
lows. We label spins by two indices; in addition to cluster 
index i, we have index within a cluster a: Sj Q . Similarly, 
interactions are denoted by Ji a ,k/3 which can be either 
ferromagnetic or antiferromagnetic. Restrictions on signs 
of interactions will be explained later, while as before we 
require | Ji a ,kp\ > J- Within each cluster, the graph of 
interactions is a connected tree. Between clusters there 
is at most one link, and the link between clusters i and k 
is present if and only if there is an edge joining vertices 
i and k. This construction is illustrated in Fig. [21 for a 
simple graph. 
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FIG. 2: An example graph with 4 vertices (left) is mapped 
onto another graph with vertices lying on the grid (right). 
Each vertex is replaced by clusters on the right. "+" and 
"— " signify values of corresponding Ti a ; solid and dashed 
lines represent ferromagnetic and antiferromagnetic interac- 
tions respectively. 

The fact that the graph of interactions within a clus- 
ter is a tree means that when considered separately, each 
cluster has a doubly degenerate ground state. We iden- 
tify them with values of fictitious coarse-grained spin Si. 
Value Si = —1 shall indicate that vertex i is included 
in the independent set. The individual spins take values 
Si = TiSi for appropriate 7j. Clearly 
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(6) 



In principle, values of Ti a can be deduced from Ji a j.p 
via the last equation, but in practice this will not be 
necessary since their values will be evident from the con- 
struction. 

Between different clusters an interaction must be ef- 
fectively antiferromagnetic. If there is a link between Sj Q 
and Skfj, we require that 

Jia,kpTi, a Tk,l3 < 0, (7) 

so that Si and Sk are coupled antiferromagnetically. 
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Consider the Hamiltonian 

Hq = — ^ ^ Jia,i(3Si a Sii3 — ^ Jia,kpSi a Sk0 
i (ia.i/3) (ia,kfi) 

^ faitySia- (8) 
ia 

Here in the first term the double summation is over the 
clusters i and over all edges (ia, i(3) within each cluster. 
In the second term the summation is over all pairs of 
distinct clusters i ^ k connected by an edge (ia,k(3). 
Individual fields h ia in JSJ can be adjusted so that Hq is 
rewritten (up to a constant) as Hq = H± + Hi with 

Hi = ^ ^ \Jia,ip\0- — TiotTipSiaSifi), (9) 
i (ia.ip) 

H 2 = |^ia,fc/3|(l - r ia S ia )(l - TkpSkp). (10) 

(ia,k(3) 

(it can be verified that all signs above are correct given 
the inequalities © and (7J|). In the above equations both 
Hi > and H2 > 0. The ground state of Ho corresponds 
to Hi = and H 2 = 0. The first equality Hi = 
implies Ti a Si a = TipSip for all a, (3 or, equivalently, Si a — 
Ti a Si for all a. With this constraint in mind, H 2 can be 
rewritten into a familiar form J2ia k/3 \ Jia,kp\(l ~ Si)(l — 
Sk)- Then H 2 = implies that no two vertices i, k that 
are joined by an edge belong to the independent set at 
the same time (Si = Sk — — 1 is forbidden). 

We now estimate the gap between ground states and 
excited states, li Hi ^ it must be at least 2\Ji a ,i/3\', 
and if H 2 ^ it is at least ^\Ji a ,kp\- Therefore the gap 
is at least 2 J. The degenerate ground states of Hq corre- 
spond to all possible independent sets. This degeneracy 
is (partially) lifted by adding a term to Hq 

Hq = Hi + H 2 + V, l/=^^ Sl0 r l0 , (11) 

i 

where is some arbitrary index within a cluster. This 
corresponds to adjusting h' i0 — hiQ — ^TiQ. Alternatively, 
rather than pick up a particular spin ia we can adjust 
h' ia — h' ia — Ti a , where fii is the size of cluster i. Then 

V= JJ2 S -^- (12) 

ia 

In either case, the net effect is that the extra term takes 
a form ^ Si, which favors independent sets with the 
largest cardinality (number of vertices with Si = — 1 is 
maximized). We need only to verify that the perturba- 
tion V is small enough to not mix the ground states and 
excited states of Hq. If cither Hi or H 2 is positive, then, 
for all spins within the corresponding cluster i we set 
Sia = T ia (this is equivalent to setting Si = +1). This 
decreases the energy by at least 2 J as we have already 
estimated. On the other hand, the energy increase due 



to degeneracy-breaking term is at most J. We conclude 
that the minima of Hq involve only independent sets, and 
that the energy gap is at least J. 

Once we have constructed a classical Hamiltonian Hq 
(|9"|)- (|12fl the ground state of which corresponds to the 
solution of NP-complete MIS problem, we formulate a 
quantum Hamiltonian H Q for use in the AQC by re- 
placing in Hq all classical spins Si a with operators <5f Q 
(cf. (J2J) and adding a driver term — r^ i=1 af . 



III. EMBEDDING PLANAR GRAPH IN A 
LATTICE 

As explained in the previous section, only planar 
graphs have to be considered. Also note that only con- 
nected graphs can be considered; for disconnected graphs 
we may solve the problem for each connected component. 
Every vertex of the planar graph is mapped onto a cluster 
of vertices and edges forming a tree. This allows vertices 
far from each other to be joined. Also we must always 
preserve an antiferromagnetic character of interactions 
between the vertices to satisfy the inequality J7J) . 

An arbitrary planar graph can be mapped onto an 
appropriately sized regular lattice of spins with near- 
est neighbor interactions, as long as such system has 
a frustration (loops with odd number of antiferromag- 
netic couplings must exist) . This is a necessary condition 
since the ground state of a non-frustrated system is triv- 
ially obtained, whereas a ground state of NP-complete 
problems typically has an exponentially large degener- 
acy, with degenerate states unrelated to each other by 
symmetry transformations. 

For illustrative purposes we now embed a graph in 
a triangular lattice with antiferromagnetic interactions 
(so that each plaquette is frustrated), but such that ev- 
ery antiferromagnetic coupling can be switched off by 
an external control. A straightforward construction of a 
graph embedding is via iterated procedure. Let's label 
the vertices with numbers 1 ... N subject to following 
constraints: a) vertex i is necessarily adjacent to at least 
one vertex j such that j < i, and b) if a set of vertices 
ji ,j 2 ,... j n , i forms a face of the planar graph and all 
jk < i, all vertices that appear inside that face must 
have numbers smaller than i. Once we constructed an 
embedding of subgraph H of graph G, we augment H 
with a new vertex i connected to it. We extend the clus- 
ters corresponding to the vertices adjacent to i and link 
them to i. A completely blind procedure may not always 
work. Let us call the vertices from H "eligible" if they 
are adjacent to the vertices from G that do not belong 
to H (G\H). Links from eligible vertices to some new 
vertex will have to be drawn at later stages. 

In the case of a non-convex embedding of a subgraph H 
we may encounter a situation where there are not enough 
qubits in 2D lattice to draw all links from the new vertex 
to eligible vertices (see Fig|3J). We can guarantee that 
this never happens if we use the following guidelines: we 
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always expand clusters corresponding to eligible vertices 
so that the embedding remains convex, eligible vertices 
lie on the perimeter of this convex region, and a minimum 
distance between eligible vertices along the perimeter is 
maintained. For the triangular lattice we make sure that 




enough space) 



FIG. 3: An example of the case where a graph's embedding 
has a non-convex shape and some eligible vertices (black cir- 
cles) lie inside of the cavities. In this case there may not 
be enough space to connect all eligible vertices to outlying 
vertices. A gray area represents a part of the graph that is 
already embedded. 

the convex object always has the form of a triangle. The 
distance between the eligible vertices along the perimeter 
must be at least two. Adding a vertex that is joined to 
subgraph H by only one edge increases the side of the 
triangle by 1 (see Fig.QJ, and adding a vertex joined by 
two or more edges to subgraph H can increase the side 
of the triangle by 2 for the worst-case scenario depicted 
in Fig. 




FIG. 4: All vertices on the top side are coupled ferromag- 
netically to the new side. A new vertex is introduced that is 
coupled antiferromagnetically [29|. 

We now explain Figs. 01 and in more details. Large 
black circles lying on a boundary of the triangle describe 
the eligible vertices; connections to them may have to 
be made at a later stage. In Fig. ^ one vertex is added 
joined to the current subgraphs by one edge. It must be 
coupled antiferromagnetically to a cluster corresponding 



new vertex 
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FIG. 5: A new vertex is connected via antiferromagnetic links. 
All other vertices coupled ferromagnetically to new sides |29ll . 



to neighboring vertex 29] . 

A "copy" of the neighboring vertex is made since more 
vertices may connect to it later on. The cluster associated 
with that vertex is expanded to include that copy. Simi- 
larly copies of all eligible vertices on the same side of the 
triangles are added to corresponding clusters. Note that 
copies of vertices are coupled ferromagnetically, whereas 
new vertex is coupled antiferromagnetically This guar- 
antees that for all large black circles the corresponding 
values of Ti a = 1 and the inequalities ©,0 are satis- 
fied. In Fig. [S] a new vertex is connected to 3 vertices 
of the subgraph H . Copies of left and right vertices are 
made by ferromagnetic couplings. No copy of the middle 
vertex is made - it cannot be connected to new vertices 
as long as we add vertices in correct order. Fortunately, 
under the iteration in both cases a minimum distance of 
2 between eligible vertices is maintained. 

Let us now make some estimates. For a connected 
planar graph with N vertices and M edges the following 
inequality holds: 

N-l < M<3(N-2). (13) 

Due to this tight bound, N alone is a good measure of 
the size of a graph. Let us now estimate the number 
of qubits needed in our construction. Adding each con- 
secutive vertex may increase the side of the triangle by 
2. A single vertex is represented by a single qubit. A 
triangle of size 2(N - 1) has N(2N - 1) < 2N 2 qubits. 
This presents an upper bound on the number of qubits 
needed to model a planar graph with N vertices. Note 
that we have used qubits quite liberally. We hope that 
far smaller number of qubits may be required for prac- 
tical implementation, possibly a number scaling linearly 
with N. Minimization of this number is an interesting 
topic for further study. 
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IV. ADJUSTABLE COUPLINGS VS. FIXED 
COUPLINGS 

The construction that we provided requires the pro- 
grammability of nearest-neighbor interactions (i.e., the 
ability to switch them off). Embedding of NP-complete 
problems can only be done into the lattices that pos- 
sess frustration. An antiferromagnet on a triangular lat- 
tice is always frustrated and therefore only the ability to 
switch off antiferromagnetic interactions is necessary in 
this case. However for a square lattice, switching between 
ferromagnetic, antiferromagnetic and zero couplings is 
necessary. This is because for a square lattice neither 
ferromagnet nor antiferromagnet has a frustration that 
requires the presence of interactions of both signs. 

Antiferromagnetic triangular lattice with adjustable 
interactions can be implemented in superconducting QC 
by one of devices described in 0] and (2^]. The idea to 
use them for implementing AQC was first put forward 
by Kaminsky, Lloyd and Orlando |2(j. However the ap- 
proach described in cannot be used when switch- 
ing between 3 states (ferromagnetic, antiferromagnetic 
and zero coupling) is needed. Moreover, the switching 
is not precisely between antiferromagnetic and zero cou- 
plings but rather between strong antiferromagnetic and 
weak antiferromagnetic couplings. While this is adequate 
for small graphs, for graphs that are sufficiently large 
(O(j on / Jos)), presence of weak antiferromagnetic cou- 
pling can lead to appearance of spurious minima unre- 
lated to the correct solution of MIS problem. This ef- 
fect can be compensated by grouping vertices in larger 
clusters and grouping links between qubits in these clus- 
ters in larger bundles. This approach is analogous to 
using thick "wires", which increases J n/JoS ratio for ef- 
fective interactions. This approach, however, is ridden 
with complexities. And most importantly, as was argued 
in Introduction, in many physical systems it is crucial to 
have an AQC architecture that does not rely at all on 
tunable interactions. This is precisely what we propose 
to accomplish in this paper. 



A. Fixed interactions 

The general idea of our approach is the following. In- 
stead of "deleting" couplings between nodes of the reg- 
ular lattice we propose to delete vertices. By applying 
a large (effective) magnetic field to a particular qubit i 
we can polarize it along the z direction (to value +1). 
Magnetic fields for its neighbors are adjusted as follows 



hi 



(14) 



The net effect is the same as if qubit i was completely ab- 
sent (not coupled to its neighbors). Since this operation 
can be applied to any qubit, we can start from a regular 
lattice and by selecting a set of qubits to be deleted ap- 
propriately, we "carve out" a desired circuit. Note that, 
although we so far assumed that spins at deleted vertices 



are completely polarized, corresponding to infinitely large 
magnetic fields, in practice we only require that spins be 
completely polarized in the ground state of the classical 
Hamiltonian. For that purpose it is sufficient to use mag- 
netic fields on the order of J, defined as J = min/^&\ J^fc. 

We identify several solutions for implementing this. 
First general approach is to use one lattice where individ- 
ual nodes are "deleted" to simulate another lattice with 
switchable interactions. The simulated lattice will have 
a larger lattice constant. Therefore a larger number of 
qubits is required. As a first example of this approach we 
take a square lattice with a particular periodic structure 
as depicted in Fig. El 111 this figure solid lines correspond 
to ferromagnetic and dashed line to antiferromagnetic 
couplings. 



FIG. 6: A patterned lattice in which 50% of elementary square 
loops are frustrated 

In Fig. we divide all qubits into four sets. "Deleted" 
qubits are represented by empty circles; the severed cou- 
plings are represented by dotted lines. "Working" qubits, 
depicted by large black circles correspond to the nodes 
of triangular sublattice. Unused "auxiliary" qubits rep- 
resented by small black circles effectively "pass on" in- 
teraction to their neighbors. The "control" qubits repre- 
sented by gray circles are used to selectively enable or dis- 
able antiferromagnetic couplings between working qubits. 
Deleting gray control qubit switches off an interaction be- 
tween two working qubits that otherwise are effectively 
coupled antiferromagnetically. As a next step we apply 
our construction for embedding a planar graph into a tri- 
angular lattice. Note that since for every working qubit 
we use 2 auxiliary, 3 control qubits and 2 deleted qubits, 
our estimate for the number of qubits is correspondingly 
increased by a factor of 8. 

For illustrative purposes we give another example 
shown in Fig. [S] This time a square lattice in which ev- 
ery plaquette is frustrated is mapped onto another square 
lattice with a lattice constant that is 3\/2 times larger. 
Now the control qubits can be used not only to switch 
off the interaction but also to choose its sign. Couplings 
between a pair neighboring working qubits is controlled 
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FIG. 7: Mapping onto triangular lattice with antiferromag- 
netic interactions. Interactions are switched off by polarizing 
gray qubits. 



by two control qubits. Deleting both control vertices 
severs the interaction between the two working qubits, 
while deleting only one makes an interaction effectively 
ferromagnetic or antiferromagnetic, depending on which 
vertex was deleted. 
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FIG. 8: Mapping onto a square lattice in which nearest neigh- 
bor interactions can be made ferromagnetic, antiferromag- 
netic or can be severed. 



Yet another approach is to do graph embedding di- 
rectly on a lattice. We choose a particular square lattice 
with frustration where each plaquette has 1 edge of one 
type (ferromagnetic or antiferromagnetic) and 3 edges of 
another type. As before, we construct the embedding in- 
teractively. Assuming that we constructed an embedding 
for a subgraph of the original graph so that all eligible 
vertices lie on the side of the square, we couple the new 
vertex antiferromagnetically to some of these eligible ver- 
tices. This can always be done since we are free to reroute 
the connections at will. An example is shown in Fig. [5] 
Control qubits as before are represented by gray circles. 
For left and right edges we delete either gray qubit. For 



the center edge we either delete qubit 2 or both qubits 
1 and 3. A simple counting argument shows that since 
there is a total of 8 choices for control qubits, exactly 
one will correspond to all couplings being effectively an- 
tiferromagnetic. This procedure increases the side of the 
square by 6 in the worst case, taking the total number 
of qubits necessary to 36A^ 2 in the worst case. As be- 
fore we stress that this is the worst-case scenario and the 
estimate will be lowered if better layout algorithms are 
used. 
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FIG. 9: In the center, either top gray qubit or two bottom 
gray qubits are deleted. On the sides either of two gray qubits 
is deleted. Since each elementary square is frustrated it is 
always possible to make all links antiferromagnetic. A gray 
area represents a part of the graph that is already embedded. 



V. FAULT-TOLERANT COMPUTING 

For nanoscale computing architectures it is important 
to take into account irregularities of the crystal lattice. 
For classical processing units in molecular computing the 
following solution was proposed A massively paral- 
lel array of simple computational devices is used, but in 
contrast to standard architectures, interconnects between 
devices are redundant (note that both devices and inter- 
connects can be defective). The presence of redundant 
connections allows to reroute the system once defective 
elements are identified. Perfectly working system was 
implemented despite the fact that 10% of elements were 
defective. 

The present AQC architecture, in which individual 
qubits arc "deleted" by polarizing them is perfectly 
suited to address the problem of fault-tolerant compu- 
tation. In our model only a small subset of nearest- 
neighbor couplings are utilized. From this standpoint, 
the system has large redundancy of links. A computer 
program, given database of defects, can achieve routings 
that bypass defects altogether. 

This can be easily illustrated as follows. Assume for 
simplicity that only the following defects can occur: (i) 
a particular nearest-neighbor coupling is too small in ab- 
solute value, \Ji 7 k\ < J for some threshold J (see Eq. (JSJ) 
and discussion after it); (ii) a particular nearest-neighbor 
coupling has a wrong sign, ferromagnetic when it should 
be antiferromagnetic or vice versa. In either of those 
cases we mark one of the qubits as "defective", delete 
it and route all links around it as shown in the Fig. 1101 
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Similarly, we can route around qubits that are themselves 
defective. This can be done as long as individual defects 
are isolated and their concentration is not too big. In 
general, it should be possible to reroute around the de- 
fects as long as their concentration does not exceed the 
percolation threshold. 

defect 




FIG. 10: If a it is always possible to route the links around 
single defect as shown in this figure. A gray area represents 
a part of the graph that is already embedded. 



Moreover, it is possible to do the routing algorithm 
even when dealing with an entirely random pattern. The 
number of possible routings increases exponentially with 
their length and with high probability we can always find 
at least one routing that makes couplings effectively anti- 
ferromagnetic. We expect that this can be done with the 
number of qubits that differs only by a constant factor 
compared to the case of regular pattern. 



VI. 



IMPLEMENTATION 



There exists a number of requirement for the imple- 
mentations of the above AQC architecture for solving 
NP-complctc problems. Firstly, the underlying spin lat- 
tice must have frustration. Secondly, an interaction be- 
tween spins must be of the Ising type, oc sfSj. Thirdly, in 
our analysis above we assumed only a nearest-neighbor 
Ising coupling. The above requirements can be satisfied 
for qubits coupled via the Heisenberg exchange interac- 
tion, superconducting qubits, and, except for the last 
condition, in many systems with "always on" magnetic- 
dipole and electric-dipole interqubit coupling. 

To see this we consider the Hamiltonian H12 for mag- 
netic coupling between the two spins- 1/2 located in the 
plane xy and subject each to a strong static magnetic 
field z B z a (a = 1,2). Also each spin is resonantly driven 
by a weak ac magnetic field B^(t) = x B^ cos(ft a t) ap- 
plied in the orthogonal direction. Assuming that the dif- 
ference between the values of the Zeeman splitting for the 
two spins is much greater then their magnetic coupling 



we obtain in the rotating wave approximation 24, 32] 

2 

(15) 

h* = l(m a - 7a B z J, T = -]p a BZ, 
D 12 (d) = , J(d) cx exp(-d/a* ) , 



where 



l/i.l, |r,-| « n 3 , \D 12 (d)\, \ j(d)\ « fi|n x - n a |. (ie) 

In CSJ we neglected small terms proportional to 
which amount to corrections quadratic in Di2/h(£li — 
O2) I and J /U{Q,\ — VL2)\, The coefficients Di2(d) and J(d) 
are magnetic-dipole and Heisenberg exchange constants, 
respectively. Both of them are positive and correspond 
to antifcrromagnetic interaction between spins. There- 
fore if spins are located at the vertices of a triangular 
lattice the later will be frustrated. 

By adjusting the amplitudes of the resonant driving 
fields Bj one can set the effective transverse field Tj = T 
for all spins except for the "deleted" spins where there is 
no ac driving (Tj — 0). Also each ac field Bj C (t) is not 
a sequence of short pulses but rather a continues wave. 
During the AQC the field detuning from resonance hj 
is fixed while the field matrix element T = T(t) is de- 
creased adiabatically slow in time. This can be achieved 
using various methods for a single spin qubit control, such 
as modification of the electron g-factor [33], and others 
that are discussed in the context of the spin-based QC in 
quantum dots [2^, |^, 0| and shallow donors [24| . 

Consider now a triangular lattice of charge qubits 
where each qubit a is encoded by the two orbital states 
with the wave functions ^„(x, y, z) (n = 0,1) and the 
difference between the energy levels equals AEj . We as- 
sume that the only nonzero dipole matrix elements for 
each qubit correspond to a z-component of the electric 
dipole vector, p J nm = — e(^ J n \z\^ J m ) . Consider now a pair 
of qubits 1 and 2 located in xy plane on a distance d from 
each other and assume that the strength of the electro- 
dipole interaction between them D 12 (d) -C AE\ — AE 2 . 
Then the truncated 2-qubit interaction Hamiltonian has 
the form H.12 = D\2(d)s\s2, where s| are z-components 
of effective spin- 1/2 operators acting on qubit states and 
Di 2 (d) = (poo — Pn) 2 /(2 hd 3 ) > is a constant of an 
electro-dipole interaction (cf. [l3)- This interaction is of 
an "antiferromagnetic" type and will lead to frustration 
on a triangular lattice of qubits located in xy plane. 

In a usual QC setting we assume that a control elec- 
trode is located near each qubit and apply a time- 
modulation of the potential of the local electrodes that 
produces weak microwave fields with the frequencies 
Qj ^> \Qj — AEj\/h. The z-component of each field 
will drive resonantly the transition between the corre- 
sponding qubit states in a near-field regime. Then in the 
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rotating wave approximation H.12 will take a form i|15|) 
(J = 0) with hj and T being, respectively, a detuning 
from the resonance and a non-diagonal matrix element 
of the field for a qubit j. By fixing {hj} and slowly 
reducing the value of T to zero one can implement the 
AQC architecture similarly to the previous case. The 
example given above is relevant for the charge qubits en- 
coded by electron states lithographically confined on a 
surface of liquid helium , and also by laterally quan- 
tized electron states in sin gle- electron quantum wells in 
GaAs heterostructures [l9ll21| . 

A main limitation of the proposed architecture for sys- 
tems with dipole-dipole interactions is that it takes into 
account only nearest-neighbor interqubit couplings. We 
note however that in our approach each vertex of the 
graph in MIS problem is embedded in a cluster with a 
large number of qubits. This picture may effectively sup- 
port the nearest-neighbor coupling approximation. The 
detailed analysis of this problem will be done elsewhere. 

Another challenge of the implementation schemes con- 
sidered above is that frequencies of resonant driving Qj 
must be different for different qubits and therefore each 
qubit has to be driven with its own ac field (in near-field 
regime). For frequency range < 10 12 Hz this can be done 
by modulating the gate bias near each qubit. This con- 
dition can be satisfied for electrons on helium 0], f° r 
electronic states in broad quantum dots |2J and for spin 
qubit control via the g- factor modulation |33t l3^ | . 

However for qubits encoded by orbital states of a shal- 
low donors in GaAs [35j with the long decay time (~ 
350ns) the intra-qubit frequencies AEj/h are in the THz 
range and local ac driving is highly problematic. One 
could consider in this case a global driving with the THz 
field containing a range of frequencies within the spectral 
window that overlaps with the intra-qubits frequencies 
{AEj/h} so that each qubit will be driven resonantly by 
its own frequency component of the field. Fixed detun- 
ing hj can be produced locally for each qubit using a 
quadratic Stark effect from the field of the control elec- 
trode. 



VII. CONCLUSION AND OUTLOOK OF 
FUTURE WORK 

In the proposed AQC architecture the interactions be- 
tween qubits are never turned off and never re-adjusted 
for solving any given instance of an NP-complete prob- 
lem. Moreover, the problem instances can be embed- 
ded into the underlying Ising lattice even in the pres- 
ence of site imperfections and random distribution of the 
magnitudes of the Ising coupling coefficients. Also in 
the proposed scheme there is no need to perform NMR- 
type refocusing sequences of fast qubit gates that are 
used in conventional quantum computing (QC) architec- 
tures |2J, |36( and must satisfy the stringent conditions 
reflecting the multiple time scales in the material sys- 
tem. In contrary, in our approach the computation pro- 



cess is guided with the slow continues-time variation of 
the parameters of single-qubit Hamiltonians avoiding un- 
wanted resonances with bulk excitations. These consid- 
erations make our analog scalable AQC architecture es- 
pecially advantageous for solid state QC implementations 
that use "alwa ys o n" magnetic- [23l . l24l.l25f and electric- 
dipole [H E lH HJ El E3 interactions for entangle- 
ment generation. This architecture is also of interest for 
QC implementations that use the Heisenberg exchange 
interaction, such as electron spin qubits in quantum dots 
[25l l3l|. It will allow to avoid an inter-qubit "J" gate, 
or any electrical control over the wavefunction overlap, 
hence making a gate lithography much simpler and re- 
ducing the sensitivity to electrode noise. 

For future work we leave the mapping of an instance 
of the Maximum Independent Set problem onto an ar- 
bitrary planar graph, edges of which represent the cou- 
plings between qubits. It would allow us to work with 
completely random systems that do not have any peri- 
odic structure. Though the problem of deciding whether 
a certain instance can be mapped onto an arbitrary pla- 
nar graph is likely NP-complete, we only need to find a 
good approximation algorithm, which can be designed to 
take polynomial time and guarantee to find such a map- 
ping for sufficiently large random planar graph with high 
probability. 

We chose to implement a Maximum Independent Set 
problem because the problem uniquely permits fluctua- 
tions in the magnitude of J^fc, as long as the sign is un- 
changed. For practical problems involving the solution of 
Constraint Satisfaction problem the problem is often first 
transformed into an instance of a Maximum Independent 
Set. Although a mapping from Constraint Satisfaction 
to Maximum Independent Set exists, it may be benefi- 
cial to work with a Constraint Satisfaction problem from 
the outset. We have described earlier that every clause 
that involves three variables can be replaced by the gad- 
get in the form of planar graph. Mapping this directly 
onto Ising model will reduce the number of extra vertices 
and edges necessary. Maximum Independent Set problem 
limits our possibilities as it assigns equal weight to ev- 
ery vertex that appears in the independent set. Relaxing 
this constraint leads for more efficient architectures for 
simulating a clause or for a gadget that replaces crossing 
edges. 

Recently Oliviera and Terhal developed a mapping 
of an arbitrary quantum circuit onto a square lattice 
with nearest neighbor interactions of three possible types: 
sfsf , s\§2, sl§2, or no coupling at all |32|- This opens 
door for emulating an arbitrary quantum computing al- 
gorithm as AQC with all spins on a 2D square lattice. 

The construction done in [37j is rigid: every quantum 
circuit corresponds to a particular lattice. It is interesting 
to explore if our approach of polarizing qubits to "carve 
out" circuits out of standard pattern can be applied in 
the case where couplings sf sf and sfsj appear alongside 
sf s|, and to find out what might be the optimal standard 
pattern. In addition it is interesting to explore if the same 
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can be accomplished starting off from a completely ran- 
dom pattern. If a standard pattern is sufficient (provided 
we can adjust fields acting on individual qubits), it may 
be a very practical approach for solid-state implemen- 
tations of a universal quantum computer and could be 
programmed to implement Shor's factoring algorithm. 
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